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Abstract 


This work shows details of an evaluation of an observational system comprising a complementary metal-oxide- 
semiconductor detector, 60cm telescope and filter complement. The system’s photometric precision and 
differential photometric precision, and extinction coefficients were assessed through observations of Supersky flat 
fields, open clusters, standard stars and exoplanets. Photometry was precision achieved at the 0.02 mag level, with 
differential photometry of 0.004 mag precision. Extinction was found to agree with previous studies conducted at 
Xinglong Observatory. Ultimately, the results demonstrate this observing system is capable of precision scientific 
observations with a charge-coupled device across the optical wavelengths. 
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1. Introduction 


Since the application of charge-coupled devices (CCDs; 
Boyle & Smith 1993) to astronomical observation, there has 
been a significant increase in observation methods. Later 
with developments in modern semiconductor technology, 
complementary metal-oxide-semiconductor (CMOS) sensors 
appeared. However, early CMOS technology was not suitable 
for scientific astronomical observation due to the high readout 
noise, and the lower dynamic range compared to CCD 
technology (Gallaway 2016). Furthermore, the independent 
readout of each pixel can potentially lead to non-uniformity in 
the readout levels between pixels. Thus, early CMOS 
technology was primarily used for non-scientific astronomical 
imaging. 

The advent of scientific CMOS image sensors led to the 
initial application of CMOS cameras in amateur astronomical 
observation (Fossum 1997; Bonanno et al. 2003). Advance- 
ments in CMOS manufacturing processes have been more and 
more significant. Currently, only a small difference exists 
between the readout noise of CMOS devices and CCDs. 
Further, the dynamic range of CMOS sensors has expanded 
considerably. In addition, the emergence of back-illuminated 
CMOS chips has substantially increased quantum efficiency 
(Bigas et al. 2006; Li et al. 2006). The readout noise of CMOS 
devices approaches that of CCDs. As a result, the difference 
between the performance of CCD and CMOS image sensors 
continues to diminish. 

Although CMOS applications are popularly used by amateur 
astronomers, CCDs are still widely used in professional 


astronomical observations. However, CMOS has been 
employed in some sky survey projects, such as the Argus 
Optical Array (Law et al. 2022), the Large Array Survey 
Telescope (Ofek et al. 2023), or the next-generation telescope 
of the ATLAS project (Tonry et al. 2018), which will be 
installed at the ATLAS-Teide observatory on Teide mountain 
(Licandro et al. 2023). Some of them use CMOS because of its 
low cost and ability of wide field-of-view and high signal-to- 
noise ratio (SNR) measurements. In fact, there are also 
differences in performance between different brands of CMOS. 
A considerable amount of testing has been done on CMOS 
(Karpov et al. 2020; Qiu et al. 2021; Alarcon et al. 2023), 
which has revealed the following unique problems such as the 
salt and pepper noise that are not seen with CCDs. However, 
these tests are solely on performance, and the astronomical tests 
of CMOS, such as photometric accuracy and extinction 
coefficient, are rare. 

In this article, we test the CMOS product SONY IMX455 
chip, which is back-illuminated. In Section 2, the observation 
system is introduced. In Section 3, the processing of the data 
reduction is described. In Section 4, the results and analysis of 
astronomical observations are introduced. Section 5 presents 
the summary and conclusions. 


2. Observations 
2.1. Observation System 


The observation system used in this article includes a 
telescope, CMOS camera and filters. 
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Table 1 

Parameters of the CMOS Camera 
Features Specifications 
Sensor SONY IMX455 
Pixel number 9576 x 6833 
Max full frame rate 2 fps 
Pixel size 3.76 um x 3.76 um 
Imaging area 36 mm x 24 mm 
A/D conversion 16 bit 
Full well 51400e— 
Readout noise 1.2 ~ 3.5e7 


Dark current noise 0.003e~ s™' pix! at 0°C 

Shutter Type rolling shutter 

Fill factor 100% 

Operating temperature Below ambient temperature 30°C ~ 35°C 


The telescope used in this observation is the 60 cm telescope 
at the Xinglong Observatory of National Astronomical 
Observatories, Chinese Academy of Sciences (NAOC). The 
telescope has the following parameters from the official 
observatory website: 


1. Aperture: 60 cm 

2. Field of view: 18’ x 18’ 

3. Focal ratio of primary mirror: F'/3.28 

4. Corrected focal ratio of primary focus: F/4.23 


The 60cm telescope achieved first light in 1964 as a 
technical test telescope for the 2.16m telescope at Xinglong 
Observatory. It was fully commissioned in 1968 after a series 
of upgrades. Since 1974, the telescope has been used for 
photometric observations of variable stars such as eclipsing 
binaries and pulsating stars. From 1995 to 2000 it was used 
primarily for a supernova search survey. Since 2000 it has been 
focused on variable star observations. 

The IMX455* (Alarcon et al. 2023) is a back-illuminated 
CMOS sensor, which is installed in the ZWO ASI6200MM 
PRO camera in this observation system. Table 1 shows some of 
the parameters of the camera, and detailed information can be 
referenced at the official website.” We can see that the dark 
current noise is only 0.003e s~! pix! at 0°C, which means 
that a 5 minute exposure would result in only 0.9e dark 
current noise, which is negligible. Pixel size is 3.76 um. Pixel 
scale is 0” 306 pixel ' for our testing system. 

We used the Johnson-Bessell filter system (Bessell 1990). 
Only BVRI are applied, and their corresponding central 
wavelengths are shown in the Table 2.ghh 


3 Reference Site: 
28.html. 

Reference Site: https: //www.sony-semicon.com/files/62/pdf/p-13_ 
IMX455AQK_BQK_ ALK Flyer.pdf. 

Reference Site: https: / /i.cmoscool.com/zwo-website /manuals / 
ASI6200MM%20Pro_Manual_CN_V1.0.pdf. 


http: //www.xinglong-naoc.cn/html/gceyq/60/detail- 
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Table 2 
Center Wavelengths Corresponding to BVRI in Johnson-Bessell filter system 
A-(am) B V R I 
Johnson 440 550 640 806(Bessell) 


2.2. Observation 
2.2.1. Observation Plan for Tests 


In this article, the tests mainly include: Supersky flat field, 
accuracy of photometry, accuracy of differential photometry 
and calculation of extinction coefficients during CMOS 
observations. 


2.2.1.1. Supersky Flat Field 


The primary goal of capturing a Supersky flat field is to 
compare it to dusk and dawn flat fields. This comparison allows 
one to assess differences between dusk and dawn sky 
illumination patterns. A Supersky flat field is obtained by 
observing low stellar density fields. A fixed exposure time was 
adopted, with the telescope being offset slightly between 
individual exposures, yielding approximately 300 images 
acquired across three bands. 


2.2.1.2. Open Clusters 


To measure the photometric accuracy, we observed four 
open clusters, NGC 6913, NGC 7243, NGC 6811 and NGC 
744, with consideration for time constraints on individual 
exposures. Open clusters were selected based on having an 
optimal density of stars within the field-of-view to balance 
sufficient numbers with avoiding excessive crowding. Addi- 
tionally, open clusters contain a large population of variable 
stars and binary systems, as noted by Abazajian et al. (2009). 
These unique features of open clusters informed our choice to 
utilize them in our calculations. Table 3 lists the specific targets 
of the photographed open clusters in our observing program. 
Open clusters serve as ideal targets to test photometric accuracy 
through calibration to literature values. 


2.2.1.3. Standard Stars 


The extinction coefficients were computed from observa- 
tions of standard stars. The method involves monitoring 
standard stars throughout the whole night from rise to set. 
Photometric standard stars have precisely determined flux 
measurements across various photometric systems (Oke & 
Gunn 1983; Landolt 2007a, 2013). Through measurements 
using a CCD camera or photometer, the brightness or flux of 
another object can be determined by comparison to the 
standards. By tracking standard star magnitudes versus airmass, 
the extinction coefficient was derived. In this work, standard 
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Table 3 
Information on the Observed Objects 


Name R.A. (J2000) Decl. (J2000) Vmag Deep 
NGC 6913 01 50 41.71 +21 45 35.89 6.6 

NGC 7243 22 15 9.12 +49 49 48.00 6.4 

NGC 6811 19 37 21.60 +46 22 40.80 7.74 

NGC 744 01 58 36.480 +55 28 22.80 8.38 

SA 20-43 00 45 42.45 +45 35 15.43 10.4 

BD+39 3312 18 04 28 +39 40 38 9.78 te 
HAT-P-32 b 02 04 10.24 +46 41 16.8 11.29 0.0244 
WASP-33 b 02 26 51.08 +37 33 02.5 8.3 0.0151 


stars were from Landolt’s catalog (Landolt 2007b, 2013). 
Table 3 lists the standard stars selected for our observations. 


2.2.1.4. Exoplanet 


To measure the differential photometric precision, we 
observe two exoplanets. Table 3 lists the two observed targets: 
HAT-P-32 b and WASP-33 b. The approach involved 
capturing single-band V exposure frames during the transit. 
Both targets were monitored continuously from one hour 
before until one hour after the transit event of the exoplanet. 


2.2.2. Observation 


Due to the limit imposed by mechanical tracking for the 
60 cm telescope, the maximum observation time for a single 
target is 30s. For the camera settings, the GAIN parameter is 
set to 100 which refers to the gain value of 0.25e /ADU, and 
the readout noise value is 1.5e . The camera was cooled 
to —10°C. 

Our observation log is listed in Table 4. The relevant 
information includes: Primary Observation Targets, Observa- 
tional time, Band, Exposure and Frame. Our main observation 
period is from 2022 September 23 through September 30. Prior 
to each night, the imaging quality was assessed and the optical 
focus adjusted to optimize image quality. Each open cluster 
was observed sequentially in four filters, with 80 images per 
band. Standard stars were monitored throughout entire nights in 
all four filters. Exoplanet targets were observed in V band alone 
for the whole night. The observation strategy employed is 
different between different targets: HAT-P-32 b was in normal 
focus, while WASP-33 b was defocused. 


3. Data Reduction 
3.1. Image Trimming 


As described in Section 2, the approximate field-of-view of 
the 60 cm telescope is 18’ x 18’. However, during observations 
it was found that the full field captured by the camera 
subtended roughly 32’ x 48’. Therefore, trimming of the 
images was required to extract the relevant field-of-view for 


Duration /minutes Type Purpose 


Open cluster 

Open cluster 

Open cluster 

Open cluster 

Standard star 

e Standard star 
186.5 Exoplanet 


photometric accuracy 
photometric accuracy 
photometric accuracy 
photometric accuracy 
extinction coefficient 
extinction coefficients 


differential photometric accuracy 
163 Exoplanet differential photometric accuracy 
Table 4 


The Observation Log Table, Listing Primary Observation Targets, 
Observational Time, Band, Exposure and Frame 


Primary Obser- 
vation Targets Observational Time Band Exposure Time Frame 
(YYYY/MM/DD) (s) 
NGC 6913 2022/09/24 B 5 80 
V 5 80 
R 5 80 
I 5 80 
NGC 7243 2022/09/24 B 5 80 
V 5 80 
R 5 80 
I 5 80 
HAT-P-32 b 2022/09/26 V 30 541 
WASP-33 b 2022/09/27 V 20 548 
SA 20-43 2022/09/28 B 20 282 
V 20 282 
R 10 282 
I 20 282 
Supersky 2022/09/29 B 20 100 
flat field 
V 20 284 
R 20 100 
NGC 6811 2022/09/30 B 30 80 
V 30 80 
R 15 80 
I 15 80 
NGC 744 2022/09/30 B 25 80 
V 25 80 
R 25 80 
I 25 80 
BD+39 3312 2023/06/21 B 4 532 
V 4 532 
R 4 532 
I 16 532 


analysis. To obtain the best illuminated field from the 
telescope, only a subsection of the full CMOS frame 
corresponding to the area illuminated by the telescope optics 
was retained for photometry. 

After the trimming, the image data size is changed to 
3354 x 3354 pixel’, and its field of view is changed to 


Research in Astronomy and Astrophysics, 24:055009 (18pp), 2024 May 


Figure 1. Image data trim; the right panel is before trim, the left panel is 
after trim. 


17/05 x 17/05. The pre-trimmed and post-trimmed planar 
fields are shown in Figure 1, where the left is the pre-trimmed 
planar field, and the box in the figure is the trimmed range. The 
right panel is the trimmed flat field. All bias frames, flat fields 
and science target images were trimmed for further processing 
and analysis. 


3.2. Bias 


This section aims to evaluate the stability of CMOS bias 
frames and investigate any temporal variations before and after 
observation with CMOS. Bias images were taken at dusk and 
dawn on different nights. Examination of the biases indicates 
the background stabilizes at approximately 510 ADU. First we 
combined bias frames using median. The median value was 
retained to mitigate the effects of salt and pepper noise 
(Alarcon et al. 2023). Differences between biases were 
analyzed to highlight any temporal variations. Table 5(a) 
shows the differences between the mean of the combined dusk 
bias and the combined dawn bias each day. Table 5(b) shows 
the differences between the mean of the 24th combined dusk 
and the other day’s combined dusk biases. 

These comparisons demonstrate the bias stability over the 
whole observation nights. Thus, we conclude the CMOS 
detector operation was stable throughout the observation 
period. During pre-processing, the dusk bias frame was 
subtracted from both the raw science images and flat fields at 
each night. Overall, the characterized bias stability supports 
reuse of a single master bias constructed from dusk frames 
only, without need to consider time-variable corrections over 
the observation run. 


3.3. Flat Field 


Flat field is essential for processing observational data not 
only to correct for the unevenness of the illumination sky, but 
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Table 5 
Bias Comparison, the Difference Value is the Numerical Difference between 
the Highest Peak and the Lowest Valley in the Resulting Image Post the 
Subtraction of the Two Bias Images 


(a) Dawn Bias versus Dusk Bias, Dusk Bias Subtracting Dawn Bias 


Date (YY YYMMDD) Difference Value (ADU) 


20220924 —0.09 

20220925 —0.048 
20220926 —0.072 
20220927 —0.054 
20220928 —0.074 
20220929 —0.007 


(b) With 24 Dusk Biases, Others days’ Dusk Bias Subtracting 24 Dusk Biases 
Date (YY YYMMDD) Difference Value (ADU) 


20220925 0.02 
20220926 0.044 
20220927 0.044 
20220928 0.063 
20220929 0.045 
Table 6 

Flat Comparison Table for Different Dusk and Dawn Dates 
Date (YYYYMMDD)/Filter B(%) V(%) R(%) I (%) 
20220924 0.5 0.4 0.4 0.9 
20220925 0.5 0.4 0.3 1.2 
20220926 0.2 0.3 0.1 1.4 
20220927 0.3 0.3 0.7 1.5 
20220928 0.4 0.5 0.6 1.7 
20220929 0.9 0.9 0.9 UNOBSERVED 


Note. The table’s filter presents the comparison result of each band, which is 
the result of the difference in the peak-to-valley value after dividing the two 
smoothed images. 


also to correct differences between the amplifiers for each 
pixel. 

Table 6 presents the outcomes of calculating the ratio of the 
median combined flat field captured during dusk versus that 
obtained during dawn for each day. This ratio represents the 
variation between the maximum and minimum of the images of 
the two flat fields after normalizing, and smoothing via the 
median filter. This metric provides a measurement of the 
illumination uniformity achieved across each band flat. 
Comparing the values among the different dusk and dawn 
flats indicates any variation introduced by changing illumina- 
tion conditions at dusk versus dawn. Table 6 provides a 
comparison within each band, revealing that the median 
difference is within 1% in the B, V and R bands, and within 
2% in the I band. 

Table 7 presents results from dividing the median value of 
the dusk flat field for each night during September 24—29. 
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Table 7 Table 8 
The Rest of the Days and 24 Dusk Flat Comparison Table SExtractor Parameters 

Date (YYYYMMDD)/Filter B(%) V(%) R(%) I (%) Parameter Specifications 

20220925 0.8 0.4 0.5 0.3 DETECT_MINAREA 50 

20220926 0.6 0.3 0.5 0.3 DETECT_THRESH 1.5 

20220927 0.7 0.6 0.4 0.4 PHOT_AUTOPARAMS 25,33 

20220928 0.6 0.5 0.6 0.3 SATUR_LEVEL 65 525.0 

20220929 0.6 0.8 0.7 UNOBSERVED MAG_ZEROPOINT 25 
GAIN 0.25 

Note. The table’s filter presents the comparison result of each band, which is PIXEL_SCALE 0.306 
BACKPHOTO_TYPE LOCAL 


the result of the difference in the peak-to-valley value from dividing the two 
smoothed images. 


These values indicate the level of consistency between daily 
flats during the observation period. As expected, flat fields 
agree to within 1% in all bands. This consistency validates the 
combined master flat field better to than 1%, but at this level it 
will not significantly impact science photometry. The temporal 
stability of the flats supports their use for corrections 
throughout the observation window without introducing 
spurious variations. 

Therefore, we took the image that had the bias correction and 
dusk flat field correction to obtain the scientific images. 


3.4. Photometry 


We used Source-Extractor (SExtractor) tools for photometric 
measurements. SExtractor’s automatic aperture photometry 
routine is derived from Kron’s “first moment” algorithm 
(Kron 1980). For details, see SExtractor’s manual (Bertin & 
Arnouts 1996).° In Table 8, we present some parameters of 
SExtractor. DETECT_THRESH represents the threshold for 
star detection, for which we have opted for the default value. 
PHOT_AUTOPARAMS denotes the parameters for automatic 
aperture photometry, for which we have also chosen the default 
settings. Additionally, BACKPHOTO_TYPE specifies the 
background calculation method, with our selection being 
LOCAL to compute the flux error using local background 
estimation. The formula used by SExtractor to calculate the 
flux error is as follows 


Fluxerr = 


(1) 


Equation (1): where A is the set of pixels defining the 
photometric aperture, g; the standard deviation (std) of noise (in 
ADU) estimated from the local background, p; the measure- 
ment image pixel value subtracted from the background and g; 
the effective detector gain in e/ADU at pixel i. Note that this 


6 Reference Site: https: / /sextractor.readthedocs.io /en /latest/Introduction.html. 


error estimate provides a lower limit of the true uncertainty, as 
it only takes into account photon and detector noise. 

With aperture photometry and target extraction we can 
obtain the light curve of the target (Huang et al. 2022). 
However, before we do this, we need to correct the exposure 
time. We need to correct the exposure time t and magnitude 
Mphot Obtained from the photometry, to determine the ultimate 
instrumental magnitude 


Minst = Mphot + 2.5 log, o(t). (2) 


3.5. System Conversion 


After the exposure time correction, we carry out the 
calibration and filter system correction. We introduce the 
reference catalog—Gaia’s synthetic photometry (Gaia-SP; Gaia 
Collaboration et al. 2023) and Gaia Data Release 3 (DR3) 
(Gaia Collaboration et al. 2023). Here we applied Equation (3), 
to make corrections and, according to different situations to 
transform the formula, to simplify the process 


MB inst = MB,Gaia—sp + cB t+ kip -X 
alls kop . (MB, Gaia—s — My Gaia-s ) 
p p 
MV int = My,Gaia—sp + cy + kiv -X 
T kay . (mp, Gaia—sp i My Gaia—sp) 
(3) 
MR int = MR,Gaia—sp + CR + Kir -X 
a kar > (My Gaia- s — MR,Gaia—s ) 
P p 
MI int = MI, Gaia—sp + cr + ki -X 


z k21 i (My ,Gaia—sp = MI, Gaia—sp) 


Here mins is the instrumental magnitude, mgaiasp is the 
magnitude with reference to the Gaia-SP (Gaia Collaboration 
et al. 2023) catalog, X is airmass, k, is the extinction coefficient 
corresponding to a certain band, kọ is the color coefficient 
corresponding to the two systems and C is the zero-point 
between the two systems. 

Flux calibrations are usually calibrated using the field star 
maps that are photographed and processed. The magnitude 
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Figure 2. This image is the division of the Supersky flat field and the dusk and 
dawn flat field, and then the image is smoothed with a box size of 100 pixel. 


calibrations used here are those that convert the instrumental 
magnitude to the true magnitude, and here we use the Vega 
system. Because of the high accuracy Gaia DR3 data, we have 
selected Gaia DR3 as the reference magnitude. 


4. Results and Analysis 
4.1. Supersky Flat Field 


This comparison evaluates the illumination differences 
between the Supersky flat field and flat field from dusk and 
dawn. The comparison reveals discrepancies in the sky 
background. Figure 2 displays the image after dividing the 
Supersky flat field and the dusk and dawn flat field, followed 
by smoothing the image to a size of 100. Table 9(a) is obtained 
by calculating the difference in terms of the peak-to-valley 
value in this figure. It presents results comparing the Supersky 
flat to dawn flats, using the same methodology for analysis. As 
flat fields were only obtained in three bands, comparisons are 
shown for those bands. The differences between dawn flats on 
each night are within one percent, matching expectations. 
Similarly, Table 9(b) compares the Supersky flat to dusk flats, 
also showing percent level agreement with the dusk reference 
each night. All flat fields show sub-percent level consistency. 
The Supersky flat field is used as an independent check of the 
dusk and dawn flat field’s ability to represent the illumination. 

Table 9 demonstrates that the flat field is stable at slightly 
under 0.01 during nights. Furthermore, each band flat field and 
the Supersky flat field remain within 1% over the entire 
observation period. This is in contrast to the Supersky flat field 
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Table 9 
Supersky and Different Day Flat Comparison 
(a) Supersky and Different Dawn Flat Comparison 
Date(Y YY YMMDD)/Filter B (%) V (%) R (%) 
20220924 0.3 0.3 0.4 
20220925 0.7 0.8 0.8 
20220926 0.4 0.4 0.5 
20220927 0.4 0.3 0.4 
20220928 0.4 0.4 0.4 
20220929 0.4 0.4 0.4 
(b) Supersky and Different Dusk Flat Comparison 
Date(Y YY YMMDD)/Filter B (%) V (%) R (%) 
20220924 0.4 0.6 0.7 
20220925 0.4 0.5 0.6 
20220926 0.5 0.5 0.5 
20220927 0.4 0.4 0.6 
20220928 0.3 0.3 0.3 
20220929 0.6 0.7 0.9 


Table 10 
Limiting Magnitudes at 3g and So for Each Band 


Filter So (mag) Accuracy of Fitting Exposure Time (s) 
B 19.152 0.042 30 
V 19.205 0.043 30 
R 18.747 0.053 30 
I 17.361 0.052 30 


for the 60cm observation field. It suggests that changes in 
observing conditions like airmass or twilight illumination are 
less than 0.01, as validated by the independent Supersky flat. 


4.2. Limiting Magnitude 


The limiting magnitude is an important indicator of an 
observing system, demonstrating the observing system’s 
capability. In this system, due to the mechanical limitations 
of the telescope, we can only use a single exposure up to 30 s. 
Below are the limiting magnitudes calculated by the open 
clusters, after calibrating and calculating the limiting magni- 
tudes in the BVRI bands and at 30 and 5ø. Figure 3 displays the 
magnitude-error plots and magnitude statistics after the 
calibration. Here the Mer is obtained from Equation (1). 

From Table 10 we see that limiting magnitude variation in 
different bands corresponds to the quantum efficiency of the 
CMOS, and the limiting magnitude at 5c in the V band is 
19.205 mag, as expected. 
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(b) V-band magnitude-error plot 
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Figure 3. The 5a SNR limiting magnitude of each band at 30 s exposures. The calibrated magnitude is on the x-axis and magnitude error on the y-axis. The blue data 
points represent the observed data, the yellow curve is a fitted curve to the magnitude error and the red markers are points with an SNR of 5c. The plot shows that the 
limiting magnitude is 19.152 mag in B-band, 19.205 mag in V-band, 18.747 mag in R-band and 17.361 mag in /-band. 


4.3. Open Clusters 


Table 3 shows four open star clusters that we observed, 
namely NGC 6913, NGC 7243, NGC 6811 and NGC 744. 
Figure 4 displays their R-band trimming images. In this section, 
we will separately introduce the imaging and processing results 
of each open star cluster. 

After completing the photometry, the reference magnitude 
we used was Gaia-SP. Based on the characteristics of open 
clusters, since in the same field, the extinction between the 
individual stars can be almost equal, we regard kı x X as a 


constant, and thus Equation (3) becomes: 


MB, int = ™B,Gaia—sp + Cg + constant 

Eg kop j (MB, Gaia—sp = My ,Gaia—sp) 
MV int = My,Gaia—sp + cy + constant 

as kay Š (MB,Gaia—sp = MY Gaia—sp) 
MR,inst = MR,Gaia-sp + CR + constant 

T kar (My ,Gaia—sp i MR,Gaia—sp) 
My int = MI Gaia—sp + c; + constant 

T k21 j (my ,Gaia—sp = MI, Gaia—sp) 


(4) 
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Figure 4. R-band images of each open cluster. 


Research in Astronomy and Astrophysics, 24:055009 (18pp), 2024 May 


* orignal 
0.107 * after select 


0.08 + ë 


Mr, gaia—sp —_ mı, gaia — sp 


2.84 2.86 2.88 2.90 2.92 
MR, inst — MR, gaia - sp 


(c) 


Mu et al. 


10 m 
a Inst üa t a 
11] GAIASP S 
12- w os 
# 
134 ce 
x Aa 3 % * 
E 14-4 A á 
ai, N 
154 a, 4 
164 A 
My kia 
174 at 
—0.75 -0.50—0.25 0.00 0.25 0.50 
Mr — Mı 
(b) 
0.044 F 
< std=0.016 
e 0.02] ° : 
£ wa 
€ ie Si uë . 3 
| gA s 
a O00) p 2. 7e «4 EN 
I : ee se e 
a —0.02 4 . po J ° OF , A 
£ ° ae 
-0.044 : 
14 15 16 17 
MR, inst 
(d) 


Figure 5. NGC 6913 R-band results. (a): The NGC 6913 R-band’s magnitude-std diagram, with red crosses representing the selected data points and blue dots 
signifying the excluded data points. (b): NGC 6913 R-band’s color-magnitude diagram, where the blue triangles represent the instrumental magnitude and red squares 
correspond to the Gaia-SP magnitude. The horizontal axis is mg — my and the vertical axis is mp. (c): NGC 6913 R-band’s magnitude—color diagram, with blue dots 
representing data points and a red line delineating the fitted line. The horizontal axis of the plot represents the difference between instrumental magnitude and Gaia-SP 
magnitude, while the vertical axis means Gaia-SP’s color, denoted as mp Gaia-sp — MI,Gaia-sp, aS described in Equation (4). (d): NGC 6913 R-band’s magnitude-error 
diagram, with blue dots representing data points. The horizontal axis represents instrumental magnitude, while the vertical axis means the difference between 


instrumental magnitude and the fit magnitude according to Equation (4). 


We can see that Equation (4) has only two variables, 
C + constant and kz, from which we can obtain their respective 
corresponding parameters by simple straight line fitting. The 
color corrections are combined using Equation (4) (Cousins & 
Caldwell 2001; Calamida et al. 2018). 


4.3.1. NGC 6913 


We performed photometry on 80 frames in the same band, 
and calculated the std of the magnitude distribution for each 
star over these 80 frame, as depicted in Figure 5(a). It displays 
the error of stars before our selection, with instrumental 


magnitude on the x-axis and the std of the magnitude 
distribution for each source across the 80 frames on the y- 
axis. Based on Figure 5(a), we obtained an internal precision of 
magnitude measurement. We then used Figure 5(a) for target 
selection, removing the variable stars, saturated stars and non- 
Zero Age Main-Sequence cases. Finally, we selected member 
stars and matched the final selected sources to the Gaia-SP star 
catalog, as depicted in Figure 5(a). 

Figure 5(b) compares the color-magnitude plot of the 
instrumental magnitude with the color-magnitude plot of 
Gaia-SP. The horizontal axis is the R—/ color, and the vertical 
axis is the R-band magnitude. Figure 5(c) is from Equation (4). 
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Figure 6. Same as in Figure 5 but for NGC 7243 R-band results. (a): The R-band’s magnitude-std diagram, with red crosses representing the selected data points and 
blue dots signifying the excluded data points. (b): The R-band’s color-magnitude diagram, where the blue triangles represent the instrumental magnitude and the red 


squares are the Gaia-SP magnitude. The horizontal axis is my — my and the vertical axis is mp. (c): The R-band’s magnitude—color diagram, with blue dots representing 
data points and a red line delineating the fitted line. The horizontal axis of the plot represents the difference between instrumental magnitude and Gaia-SP magnitude, 


while the vertical axis means Gaia-SP’s color, denoted as mg Gaia-sp — MGaia-sp, aS described in Equation (4). (d): The R-band’s magnitude-error diagram, with blue 
dots representing data points. The horizontal axis represents instrumental magnitude, while the vertical axis means the difference between instrumental magnitude and 


the fitted magnitude according to Equation (4). 


The vertical axis in Figure 5(d) is the difference between the 
corrected magnitude mp g and Gaia-SP mgaia-sp. After the final 
fit, we calculate the difference between the magnitude and 
Gaia-SP color by multiplying the color coefficient with the 
color value. The R-band instrumental magnitude is the 
horizontal axis, and the accuracy obtained from the final 
photometry is reflected in the std of the label in Figure 5(d). 
Table 11 displays the coefficients k2 R, zeros cro + constant 
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Table 11 


The Results about the Observed Open Cluster, k2 R and co,z + constant 


ID 


Corresponding to the Parameters of Equation (4) 


kor Cro + constant o 
NGC 6913 0.208 2.789 0.016 
NGC 7243 0.168 2.864 0.022 
NGC 6811 0.149 2.734 0.011 
NGC 744 0.208 2.682 0.019 
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Figure 7. Same as in Figure 5 but for NGC 6811 R-band results. (a): The R-band’s magnitude-std diagram, with red crosses representing the selected data points and 
blue dots signifying the excluded data points. (b): The R-band’s color-magnitude diagram, where the blue triangles represent the instrumental magnitude and the red 
squares are the Gaia-SP magnitude. The horizontal axis is mp — my and the vertical axis is mg. (c): The R-band’s magnitude—color diagram, with blue dots representing 
data points and a red line delineating the fitted line. The horizontal axis of the plot represents the difference between instrumental magnitude and Gaia-SP magnitude, 
while the vertical axis means Gaia-SP’s color, denoted as mp Gaia-sp — MI,Gaia-sp, aS described in Equation (4). (d): The R-band’s magnitude-error diagram, with blue 
dots marking data points. The horizontal axis represents instrumental magnitude, while the vertical axis means the difference between instrumental magnitude and the 


fitted magnitude according to Equation (4). 


and photometric accuracy that were computed for each open 
cluster. From this table, we can see that its photometric 
accuracy is 0.016 mag. 


4.3.2. NGC 7243 


The processing procedure for NGC 7243 is similar to that of 
NGC 6913. Figure 6(a) shows the magnitude and std plot 
obtained with the same processing procedure for NGC 6913. 
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Finally, by calculating the std of the data in the scatter plot of 
Figure 6(c), we obtained a photometric accuracy of 0.022 mag. 


4.3.3. NGC 6811 


The processing procedure for NGC 6811 is similar to that of 
NGC 6913. Figure 7(a) shows the magnitude and std with the 
same processing procedure as NGC 6913. Figure 7 also 
represents the same calculation process as Figure 5 for NGC 
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Figure 8. Same as in Figure 5 but for NGC 744 R-band results. (a): The R-band’s magnitude-std diagram, with red crosses representing the selected data points and 
blue dots signifying the unselected data points. (b): The R-band’s color-magnitude diagram, where the blue triangles represent the instrumental magnitude and the red 
squares are the Gaia-SP magnitude. The horizontal axis is mp — my and the vertical axis is mg. (c): The R-band’s magnitude—color diagram, with blue dots representing 
data points and a red line delineating the fitted line. The horizontal axis of the plot represents the difference between instrumental magnitude and Gaia-SP magnitude, 
while the vertical axis means Gaia-SP’s color, denoted as mg Gaia-sP — MI,Gaia-sp, 2S described in Equation (4). (d): The R-band’s magnitude-error diagram, with blue 
dots representing data points. The horizontal axis represents instrumental magnitude, while the vertical axis means the difference between instrumental magnitude and 


the magnitude obtained by fitting according to Equation (4). 


6913. Finally, we calculate the std of the data in the scatter plot 
of Figure 7(c), with photometric accuracy of 0.011 mag. 


4.3.4. NGC 744 


The processing procedure for NGC 744 is similar to that of 
NGC 6913. Figure 8(a) shows the magnitude and std plot 
obtained through the same processing procedure as NGC 6913. 
Figure 8 also represents the same calculation process as 
Figure 5 for NGC 6913. Finally, by calculating the std of the 
data in the scatter plot of Figure 8(d), we obtained a 
photometric accuracy of 0.019 mag. 
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Table 11 summarizes the results for these four open star 
clusters, where the std listed in the table represents the 
photometric accuracy. The color coefficient kz is not the true 
system conversion coefficient and incorporates the main order 
contribution of the cluster. Therefore, the photometric accuracy 
is approximately 0.02 mag (Landolt 2007a; Munari et al. 2014; 
Yang et al. 2021). 


4.4. Standard Stars 


Table 3 shows that we observed two standard stars, namely 
SA 20-43 and HD 165434, in four different bands. After 
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Figure 9. Image of the four bands of BD+39 3312 taken by our observing system. The blue dots represent the observed data, while the red line delineates the fitted 
line. The x-axis represents airmass, and the y-axis represents the difference between the instrumental magnitude and the Gaia-SP reference magnitude (residual). 


Table 12 
The Results of the Four Bands on BD+39 3312 are Obtained by Iterative 
Correction 
ky Ok Co Oco 
B 0.335 0.022 0.786 0.024 
V 0.193 0.018 1.107 0.019 
R 0.090 0.013 1.122 0.014 
I 0.070 0.017 1.011 0.018 


completing the pre-processing and photometric processing, we 
need to extract the light curves of the targets in each band, and 
then match them with the catalog. After completing these steps, 
we calculate the extinction coefficients using the system 


transformation, Equation (2). At this point, the system 
transformation formula needs to be modified. Because there 
is only one target whose color is the same, for this standard 
star, the color term kz X (Mg, Gaia — Mv,Gaia) can be regarded as a 
constant, and then the formula can be rewritten as: 


kip - X + constant 


MB int = ™B,Gaia—sp CB 
My,inst = My,Gaia—sp + Cy + ki,v - X + constant 
kr + X + constant ` 


(5) 


MR,inst = MR,Gaia—sp CR 
Miint = MY,Gaia—sp + Cr + kiz + X + constant 

We can see that Equation (5) has only two variations, 
c + constant and kı, from which we obtain their respective 
corresponding parameters by linear fitting. 
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Figure 10. SA 20-43 V-band Mesiduais VS. UT. 


4.4.1. BD+39 3312 


Figure 9 shows the airmass versus magnitude of the standard 
star BD+39 3312 in each band. The blue dots represent the 
observed data, while the red line delineates the fitted line. The 
x-axis represents airmass, and the y-axis means the difference 
between the instrumental magnitude and the Gaia-SP reference 
magnitude Mresiqual. Because the observations were made in 
summer, the range of airmass values shown here is relatively 
small. Table 12 provides a summary of the fitting results in four 
bands, with k, and co corresponding to the respective 
coefficients in Equation (5). o% and o, represent the simulated 
errors. kig is 0.335+0.022, kıy is 0.193 + 0.018, kır is 
0.090 + 0.013 and kı z is 0.070 + 0.017. 


4.4.2. SA 20-43 


Figure 10 shows the variation of the standard star’s 
magnitude with respect to the observation time. The x-axis 
represents the time in UT (hr), while the y-axis represents the 
difference between the instrumental magnitude and the Gaia-SP 
magnitude (Mesiduai). The variation pattern shows a gradual 
dimming followed by a sudden brightening, continuing to 
increase in brightness until reaching the maximum. This point 
indicates the target reaching the zenith. Subsequently, the target 
undergoes a descending phase with a continuous dimming, 
followed by a slight brightening and eventually fading. This 
pattern is observed for all targets within the field of view, 
indicating the absence of any variability during the observation 
time. Figure 11 displays the airmass versus magnitude, with the 
x-axis representing airmass and the y-axis representing the 
difference between the instrumental magnitude and the Gaia-SP 
reference magnitude (Mesiqua)- The blue triangles represent the 
ascending phase, while the orange dots signify the descending 
phase. 
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Figure 11. V-band airmass vs. magnitude plot of SA 20-43. Triangles are UP 
segments, dots are DOWN segments. 


This problem has also been discussed in the literature (Yan 
et al. 2000; Zhou et al. 2001), with two different methods for its 
correction: One is to correct for the effect of atmospheric 
extinction by including a band-independent but time-varying 
zero correction term in its calibration, which was later replaced 
by a time-varying atmospheric extinction coefficient term (Yan 
et al. 2000). The other is to use a time-varying atmospheric 
extinction coefficient term to correct for the effect of 
atmospheric extinction variations (Zhou et al. 2001). In this 
work, we used the second method to set up a time-dependent 
term for the atmospheric extinction coefficient. Since we are 
only using one target, there is no the color item. To better 
correct atmospheric extinction, we add a time-varying atmo- 
spheric extinction coefficient term. For Equation (5), based on 
the corrected relationship of the individual standard stars 
above, this becomes: 


kı = kı + f (UT). 
f(UT) =a. UT +b. 


(6) 
(7) 
™MB,Gaia—sp + CB 

(kıs + ag - UT + bg) - X + constant 

My inst = My,Gaia—sp + CV 

(ki,v + ay - UT + by) - X + constant 
MR Gaia—sp + CR i 
(kır + ar : UT + bpr) - X + constant 

MI int = MI Gaia-sp + CI 

+ (kız + ar - UT + br) - X + constant 


MB int = 


(8) 


MR int = 


Here, in reality, Equation (7) is a polynomial. We start with 
the simplest case for our analysis. By fitting the data using this 
approach, we can obtain the coefficients. Subsequently, 
through a stepwise iteration, we can transform the monomials 
into a polynomial, thereby fully correcting the observed data. 
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Figure 12. The corrected data and fitting results for SA 20-43 are presented. The x-axis represents airmass, while the y-axis means the difference between the 
instrumental magnitude and the Gaia-SP reference magnitude (mM,esiduai). The blue dots represent the corrected data, while the red line delineates the linear fit obtained 


from the corrected data. 


Table 13 
The Results of the Four Bands on SA 20-43 are Obtained by Iterative 
Correction 
Filter kı Tk Co Teg 
B 0.303 0.013 2.192 0.036 
V 0.192 0.012 2.622 0.032 
R 0.137 0.009 3.245 0.024 
I 0.089 0.006 3.826 0.028 


This yields results shown in Figure 12, where the x-axis 
represents airmass and the y-axis means the difference between 
the instrumental magnitude and the Gaia-SP reference 
magnitude (M,esiduai). The blue dots represent the corrected 
data points, while the red line delineates the linear fit of the 
corrected data. The final results are summarized in Table 13. 
The coefficients kı and co in Table 13 correspond to the 
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respective coefficients in Equation (8). op and o,, represent the 
simulated errors. kı g is 0.303 + 0.013, kıy is 0.192 + 0.012, 
kı r is 0.137 + 0.009 and k; z is 0.089 + 0.006. 


4.4.3. Comparison 


Table 14 and Figure 13 compare all our measured extinction 
coefficients with those calculated for other telescopes at 
Xinglong Observatory (NAOC). 

The results were nearly identical to others measured at 
Xinglong Observatory (NAOC). 


4.5. Exoplanets 


For exoplanets, we use differential photometry, which 
compares the light curves of stars within the field of view of 
an image, eliminating the influence of instruments and weather. 
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Table 14 
Other Working Extinction Coefficients at Xinglong Observatory (NAOC) 
Article Telescope Year kip kiv kir kia 
BD+39 3312 (this article) 60 cm 2023 0.335 + 0.022 0.193 + 0.018 0.090 + 0.013 0.070 + 0.017 
SA 20-43 (this article) 60 cm 2022 0.303 + 0.013 0.192 + 0.012 0.137 + 0.009 0.089 + 0.006 
Bai et al. (2018) 85 cm 2016 0.431 + 0.029 0.282 + 0.026 0.217 + 0.019 0.156 + 0.021 
Huang et al. (2012) 80 cm 2011-2012 0.348 + 0.022 0.236 + 0.017 0.168 + 0.019 0.085 + 0.021 
Huang et al. (2012) 80 cm 2006-2007 0.307 + 0.009 0.214 + 0.008 0.161 + 0.008 0.091 + 0.008 
Huang et al. (2012) 80 cm 2004-2005 0.296 + 0.012 0.199 + 0.009 0.141 + 0.010 0.083 + 0.009 
Zhou et al. (2009) 85 cm 2008 0.330 + 0.007 0.242 + 0.005 0.195 + 0.004 0.066 + 0.003 
Shi et al. (1998) 60 cm 1998 0.31 0.22 0.14 0.10 
Shi et al. (1998) 60 cm 1995 0.35 0.20 0.18 0.16 
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Figure 13. With other working extinction coefficients, the x-axis represents 
wavelength, while the y-axis means the extinction coefficient. The data labeled 
as BATC are sourced from Zhou et al. (2001), the data labeled 2012 are 
sourced from Zhao et al. (2020) and the data labeled 2011 are sourced from 
Huang et al. (2012). 


The real used to correct the star, the reference star and the 
brightness of reference star cannot be too different from the 
target star, and the distance between the figure and the target 
cannot be too far. Additionally, the sigma of the light curve 
cannot be too large, i.e., 0.01, to ensure that it is not a variable 
star. The number of reference stars should be small, so that the 
image is the most appropriate. A comparable star, with the 
same requirements as the reference star, is also required at the 
end for final comparison. 


4.5.1. HAT-P-32 b 


For differential photometry, we can obtain its light curve. 
Figure 14 is the light curve extracted by us after differential 
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is deltaTime (hr), representing the time span from the beginning to the end of 
the observation. The y-axis of both panels represents the relative magnitude, 


My relative- 


correction. The top panel is the light curve of the host star of 
this exoplanet, and the bottom panel is the light curve of the 
reference star. The std in the label in the bottom panel stands 
for differential photometric accuracy, with scatter of 
0.0063 mag. Then, we combine the five frames into one frame 
(binned 5) to increase the SNR in Figure 15(a). Figure 15(a) 
represents our photographic acquisition of HAT-P-32 b and its 
control star, and we can see the obvious eclipse light variation 
in the target, with scatter of which is 0.0037 mag. Figure 15(b) 
is the comparison of our light curve with those of other work 
from the literature, and they completely agree with each other, 
but the basic trend is clear. The data used in Figure 15(b) are 
described as follows. The observing equipment they used in the 
articles (Hartman et al. 2011; Nortmann et al. 2016) was the 
OSIRIS instrument of the Gran Telescopio CANARIAS (GTC) 
in long-slit spectral mode. 
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Figure 15. Comparison of HATP-32b light curve and HATP-32b light curve. (a) The eclipse light curve of HAT-P-32 b is shown in Figure 14, with data binned at a 
size of 5. (b) Comparison of eclipse light curves of HAT-P-32 b. Blue data points represent observations from this observation, i.e. the data in Fig. 1, and orange data 


points represent eclipse light curves from Nortmann et al. 2016. 
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Figure 16. The upper panel shows the observed eclipse light curve of WASP- 
33 b, while the lower panel displays the variation of the comparison stars. The 
same x-axis for both panels is deltaTime (hr), representing the time span from 
the beginning to the end of the observation. The y-axis of both panels 
represents the relative magnitude, my relative- 


4.5.2. WASP-33 b 


Since WASP-33 is a 6 Scuti variable star (Herrero et al. 
2011), this means that the light curve of the exoplanet contains 
the variation of the host star. The relationships in Figures 16 
and 17(a) for WASP-33 b are the same as those in Figure 14, 
where the scatter is of 0.0066 mag, and Figure 15(a) for HAT- 
P-32 b, taken before and after binning 5. Figure 17(a) shows 
our light curve of WASP-33 b and its host star. We can even 
see its light variation of 0.02 mag clearly, while the precision of 
the former reference star is about 0.004 mag in the best 
situation. This depth of the evolutionary eclipse may include 
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the light variations of the host star. Comparing WASP-33 b 
with data obtained by other observers using CCD observations, 
our result also agrees well with Johnson et al. (2015) obtained 
from the 0.3 m telescope at Monte Carbre Observatory and the 
0.8 m telescope at Montserrat Observatory in Figure 17(b). 

We compare the observing performance of the 60cm 
telescope using the CCD. From the articles Wang et al. 
(2017) and Wang et al. (2019), we see that they have 60cm 
frames set up for exoplanets with differential photometric 
accuracy around 0.0016—0.0063 mag. The differential photo- 
metric accuracy obtained from our measurements is between 
0.003 and 0.004 mag, which is comparable or even smaller than 
the references. Therefore, we conclude CMOS is suitable for 
observing exoplanets and can replace the CCD for this 
program. 


5. Summary and Conclusions 


In this work, the tests mainly include these items: Supersky 
flat field, the accuracy of photometry, the accuracy of 
differential photometry and the extinction coefficients with 
CMOS observations mounted on the 60cm telescope at 
Xinglong Observatory (NAOC). Based on the tests, flat field 
correction for CMOS + 60cm telescope could be better than 
1% at BVRI bands, and thus there is a change of flat field every 
day during the observation period. A flat field can be replaced 
by that of another day during our observations. Therefore, we 
conclude that the camera’s performance is stable. The 
following are our conclusions: The open clusters result in a 
photometric accuracy of 0.02 mag. The calculated extinction 
coefficient of our work is compared with other coefficients 
observed at Xinglong Observatory (NAOC) using a standard 
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Figure 17. WASP-33 b light curve we collected and comparison of the different WASP-33 b light curves. (a) The eclipse light curve of WASP-33 b is shown in 
Figure 16, with data binned at a size of 5. (b)Comparison of eclipse light curves of WASP-33 b. Blue data points represent observations from this observation, i.e. the 
data in Fig. 1, and orange data points represent eclipse light curves from Johnson et al. 2015. 


star, not revealing a significant difference. The differential 
photometric accuracy for exoplanets is 0.004 mag. 

Therefore, we can conclude that the observation results of 
the observing system in this paper are in accordance with our 
expectations and can be used in scientific observations. This 
SONY IMX455 CMOS sensor satisfies our requirements and 
can replace the CCD. 

In the future, we will test the pixel stability of this camera. 
We will also test it on other devices with exposure times of 60, 
300, 600 and 1200s. Additionally, we will study the dark 
current for each corresponding exposure time. 
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